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Abstract 

We make an attempt to clarify the role of the annihilation or "penguin" mode in 
the description of the K — ► tttt decay within the Standard Model. The attention is 
concentrated on new operators in the effective AS = 1 Hamiltonian and the viola- 
tion of factorization for mesonic matrix elements of the local four-quark operators. 
We propose a regular method to evaluate the mesonic matrix elements of K — > tttt 
transitions based on studying three-point correlators via QCD sum rules using the 
chiral effective theory as an underlying low-energy model for strong interaction. 
Matrix elements of the QCD penguin operator are calculated within this approach. 
The total "penguin" contribution is found to be relatively large that improves the 
theoretical description of the AI = 1/2 rule in non-leptonic kaon decays. 
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1 Introduction. 



At present the Standard Model (SM) |1J of electroweak and strong interactions seems to 
fit successfully all known data of particle phenomenology and provides a good prototype 
of future theory of unification. Even though the SM itself is a little unsatisfactory from 
the aesthetic point of view, it plays an essential role as a convenient tool of compact 
description of all available experimental facts. However along with many achievements 
of the SM there exist several subtle points still to be understood. These points are of 
particular interest because their resolution might lead to a discovery of new features absent 
within the SM itself and could serve as an indication on new physics beyond the SM. Before 
claiming the discovery of new physics and demanding going beyond the SM, however, the 
corresponding phenomenon should be examined thoroughly to guarantee it cannot be 
explained at the level of the SM and the proper accuracy is achieved. Unfortunately this 
task cannot presently be easily accomplished. 

One example of these subtle points within the SM is the origin for considerable en- 
hancement of the AI = 1/2 parts of the amplitudes in non-leptonic kaon decays. Though 
the AI = 1/2 rule is explained qualitatively by the influence of strong interactions the de- 
tailed quantitative description is still lacking. The closely related problem of the "direct" 
CP violation in kaon decays also requires a more precise quantitative analysis. 

Recently numerous attempts have been made to improve the theoretical description 
of non-leptonic kaon decays in SM with sufficient accuracy. The efforts are exerted both 
to extend the perturbative QCD analysis @ beyond the leading order and to account 
more adequately for long-distance effects using some models of strong interactions at low 
energy such as, for example, chiral Lagrangians ||, 1/N C expansion of many color QCD 
||, |5j or lattice simulations ||. However the present results still do not reproduce the 
experimental data and the existing accuracy is not satisfactory. 

In the present paper we study some questions connected with the annihilation mode 
of the K — > mi decay which were not considered in the previous analysis. This mode is 
important in explanation of the AI = 1/2 rule because it is purely of the AI =1/2 sort 
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and its presence is the main distinction between AI = 1/2 and AI = 3/2 amplitudes at 
least within the perturbation theory. It also generates the imaginary part of the effective 
Hamiltonian and, therefore, is responsible for the non-zero value of the e' parameter 
describing the direct CP violation in kaon decays 

At the fundamental level of the SM the strangeness changing transitions with AS = 1 
occur via the VT-boson exchange between two weak charged currents. The short-distance 
analysis of the product of weak hadronic currents after removing the VT-boson and the 
heavy quarks results in an effective AS = 1 Hamiltonian of the following form [|7|, [|], [| 

H AS =i = ^WEhW + TyiMQi] + h.c. (1.1) 
V2 i=x 

Here Gf is the Fermi constant, V stands for the quark flavor mixing matrix, r = 
— VtdVts/VudVusi z %{^) an d Vi{n) are the coefficients of the Wilson expansion subtracted 
at the point fi and {Qi\i = 1, ...,6} is the basis of the local operators containing light 
quark fields (u,d,s) only with dimension six in mass units (four-quark operators). In 
eq. (1) we have omitted the contributions of electroweak penguin operators || |10[. The 
current-current operators {Qi\i = 1, 2} do not form a close set under renormalization and 
additional so-called QCD penguin operators {Qi\i = 3, ...,6} are produced by an annihi- 
lation diagram. Since the discovery of the annihilation mode it was widely discussed 
in the literature. The present analysis however cannot be considered a complete one. 
Following points require further investigation: 

1. The QCD perturbation theory analysis of the kaon non-leptonic decay can be improved 
by using more accurate effective Hamiltonian. In the standard approach only the leading 
terms in the inverse masses of heavy quarks are taken into consideration while a proper 
account of the non-leading corrections in the inverse mass of charmed quark generated by 
the penguin-type diagrams along with the usual expansion in the strong coupling constant 
a s can be important. Analogous corrections to the K° — K° mixing have been considered 
in ref. [|11[] and turn out not to be negligible though their actual magnitude can strongly 
depend on the procedure used for estimation of the corresponding matrix element. Further 
corrections to the effective Hamiltonian appear because the top quark is heavier than other 
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quarks and W-boson. This results in an incomplete GIM cancellation of the annihilation 
diagram including t-quark. Though this question was widely discussed || the appearance 
of a new operator in the effective Hamiltonian remained beyond consideration. 

2. An achievement of a high precision in theoretical estimates for the kaon decays stum- 
bles at a necessity to calculate mesonic matrix elements of local four-quark operators in 
the effective Hamiltonian (1.1). The only method of computation entirely based on first 
principles seems to be numerical simulations on the lattice though so far even this ap- 
proach has not provided us with unambiguous estimates due to some subtleties connected 
with the description of fermions on the lattice. Meanwhile several semi-phenomenological 
techniques have been developed and applied for computation of those matrix elements, 
for example, the many color expansion of QCD 0, (5|. Their precision, however, still need 
to be essentially improved. 

Recently a new regular method to evaluate the mesonic matrix elements has been 
proposed ]12[ where the effective Hamiltonian is represented in terms of the chiral effective 
theory variables and the parameters of the chiral representation are determined via QCD 
sum rules for an appropriate three-point Green's function. In Sect. 3 we briefly describe 
the method taking as an example the calculation of the QCD penguin operator matrix 
elements which has been a subject of a few controversies j|, [5], [7[ [13], PH - This technique 
will be applied in Sect. 4 to estimate the scalar gluonium contribution to the K —* nn 
decay amplitude. 

3. The renormalization group improved perturbation theory does not take into account 
the strong interaction of the soft light quarks and gluons with the virtual momentum 
smaller than the normalization point \i ~ 1 GeV. The information on this interaction 
is entirely contained in the mesonic matrix elements of the local four-quark operators. 
Well known factorization procedure for evaluation of these matrix elements |15| accounts 



only for the " factorizable" part of the interaction ||16|| . " Unfactorizable" contributions, 
for example, those corresponding to the annihilation of a quark pair from the four-quark 
operator into soft gluons are omitted when the factorization procedure is performed. The 
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calculation of these contributions and the generalization of the matrix elements beyond 
the factorization framework can be systematically done within the approach applied in 
Sect. 3 to the calculation of the penguin operator matrix element. It turns out to be 
possible to obtain an information on the important part of the factorization violating 
contribution to the K — > tttt decay amplitudes. This possibility is connected with the 
investigation of a new K — > tttt decay channel induced by annihilation of the quark 
pair from the four-quark operator into gluons with the subsequent formation of the pion 
pair by the soft gluon cloud, i.e. the decay channel with the gluons playing the role of 
the intermediate state |17]]. Being unfactorizable this decay mode does not appear as a 
correction to some leading order contribution and can be studied by its own. This feature 
makes obtained results more accurate. 

In the Sect. 4 we study a new K — > tttt decay channel with the simplest scalar color- 
less gluon configuration forming an intermediate state. We calculate both short- distance 
(perturbative) and long-distance (non-perturbative) part of the corresponding amplitude. 
The short-distance analysis is based on the results of the Sect. 2. To obtain the long- 
distance contribution the chiral Lagrangians are used as a low-energy model of strong 
interactions and corresponding chiral coupling is derived via QCD sum rules. 

2 The local effective AS = 1 Hamiltonian beyond the 
four-quark approximation. 

We begin with the calculation of the non-leading charmed quark mass correction. Expres- 
sions for Wilson coefficients Zi(fj) and yi(fi) in eq. (1.1) in the region of QCD asymptotic 
freedom are obtained by performing the operator product expansion (OPE) of two weak 
charged quark currents. At a typical energy scale of weak decays of light hadrons these 
Wilson coefficients have been with the renormalization group technique. After removing 
all heavy particles (IV-boson, t-, b-, and c-quarks) from the light sector of the theory 
formula (1.1) corresponds to the leading order in inverse masses of these particles. The 
removal of the c-quark however is not very reliable and, in general, requires a special in- 
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vestigation it is not heavy enough in comparison with the characteristic mass scale in the 
sector of light u-, d-, and s-quarks, for example, with the p-meson mass. The non-leading 
terms in the l/m c expansion can, therefore, be important and require a quantitative 
consideration. 

The effective low-energy tree level Hamiltonian for AS = 1 transitions before decou- 
pling the c-quark reads 

HX s =i = ^V ud V: s (Q u 2 - (1 - r)Ql) + h.c. (2.1) 

where Q\ = 4(s f ,7 At g z ,)(g£,7 At <i z ,), qL(R) stands for left-(right-) handed quark. Performing 
the OPE and restricting oneself to the first order terms in a s and m~ 2 one can write the 
representation for the effective Hamiltonian in the form 

H AS=1 = H^ + H®. (2.2) 

The first addendum on the right-hand side (rhs) of eq. (2.2), corresponds to the 

leading contributions in l/m c and coincides with the rhs of eq. (1.1). Second addendum 
on the rhs of eq. (2.2), is the leading order l/m c correction |L8] 



(2.3) 



where a basis {Q\ = 1, 7} |i = 1, •••,4}) of the local operators with dimension 

eight (seven) in mass units is chosen in the form 

Q% — $ L{PaG fialvG vfi + ^yG y^G ^a-FoO^L; 

(8) A A 

Qf = s L (D 2 J)d L , Qf ] = ts L (DG ull G Ufl - G v , L DG vp )d L , 

Ql^ = Sr(G ^yO ^ V G ' a {iO a /3)dL, Q?P = S r(G ^Gy^dh, 
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= is R {G va G aix a Vil )d L , = s R (J^ + PMd L . (2.4) 

Here P M = id^ + g s A^ is the momentum operator in the presence of the external gluon field 
= A®t a , t a are the standard generators of the color group SU(3), G^ v = G a ^ u t a is the 
gluon field strength tensor, = Y, q = u d siQl^o)^ an d cr^u = Jv]/2- When eq. (2.3) 
is derived we consider u- and <i-quark to be massless, keep the first order quantities in 
strange quark mass and use the equations of motion 

sP = m s s, Pd = 0, 

[P M , G^v] = iD^G^ u = -ig s J v . (2.5) 

Straightforward calculations give the following numerical values for the coefficients C\ 
to the leading order in the strong coupling constant a s 

r (8) _ 8 ( 8) _ 16 (8) _ 4 {8) _ 2 
Gl "15' ° 2 ""15' ° 3 "~5' ° 4 "15' 

_ n r (s) _ 8 r (8) _ 2 
° 5 " U ' ° 6 _ "15' 7 ""15' 

CP=~ = ^ 7) = | CP = 0. (2.6) 

Thus we have generalized the effective Hamiltonian for AS = 1 decays beyond the leading 
order in l/m c . 

To complete our treatment of the local effective Hamiltonian the annihilation of the top 
quark should be also considered. Because the top quark is heavy enough with respect to 



the jy-boson [19j these two particles should be integrated out simultaneously to obtain the 
local effective Hamiltonian. The procedure of decoupling the heavy top quark is described 
in details in ref. || but one point has remained beyond the analysis. It is connected with 
other penguin-type operator m s Q^ = m s s R g s G d^. This quark-gluon operator was 
omitted in the early papers [|7|, § because the corresponding coefficient function vanishes 
in the first order in a s . The next-to- leading calculations [^(| has verified the smallness of 
this Wilson coefficient. However the above statement holds only under the assumption 
m t <C Mw, i-e. in the leading order in m^/M^ when the complete GIM cancellation 
takes place. In the case of a heavy top quark it becomes invalid and the quark-gluon 



operator m s (5 appears in the effective Hamiltonian already in the first order in a s . 
This additional contribution has been obtained in ref. pi] and reads 

AifW = 9± VudV:sTC ^^) ms Q^^)^ 

\ Til?' 

C ifi) = 16V (F(Xc) ~ F{Xt)) ^ Xq = m|' 



o a , o 39 



F(x c ) ~ F(0) = | (2.7) 

The renormalization group factor 

^ _ f -<m b ) V (5V2ft ^K) V (5)m f a.(A*) V <5)/% 



\a s (M^)y \a s {m b ) ) \a s {m c ) J 

where 7^ = —28/3 is the anomalous dimension of the operator m s Q^ = 
11 — |n/, n/ is the number of active quarks flavors, can be easily obtained because the 
operator m s Q^> does not mix with the four-quark operators in the leading order. 

The contributions of the u- and c-quarks to the real part of the effective Hamiltonian 
cancel one another via GIM mechanism and the operator m s Q^ mainly contributes to the 
imaginary part of the effective Hamiltonian and, therefore, can be important in the anal- 
ysis of direct CP violation. However, its contribution is suppressed numerically because 
r](n) < 1 and the function F(x) changes slowly. Indeed, at the point Aq C d = 0.3 GeV, 
fx = 1 GeV, m t = 130 GeV, we have = 0.0009, while the numerical value of the 
same Wilson coefficient of the dominant penguin operator Q 6 = — 8 X^=«,d,s('SL ( ta)( , ?i? c ^) 
is y 6 = 0.102 [|. 

Thus, the complete form of the effective Hamiltonian up to the first order in m s , l/m c 
and a s with m t ~ My/ becomes available now (eqs. (2.2, 2.3, 2.7)). The electroweak 
penguin operators have to be also included into the complete expression. 

As an example of using the above Hamiltonian we consider the K — > tttt decays. For 
this end we have to extract an information about the matrix elements of the local operators 
between the mesonic states. We start with the operators m s Q- 7 ^ and m s Q^ . These 
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operators contain explicitly the strange quark mass and, therefore, in the leading order of 
the chiral expansion they correspond to the tadpole term in the chiral weak Lagrangian. 
In detail the tadpoles will be considered in Sect. 3. Now we only note that tadpoles do 
not generate any observable effect and can be neglected in the leading order of the chiral 
symmetry breaking. 

Here a remark about the role of the operator m s Q^ in determination of the parameter 
e' is necessary. Recently an enhancement of the corresponding Wilson coefficient in the 
next-to-leading order in a s has been discovered and a significance of this operator in the 



analysis of the direct CP violation was announced However in ref. [23| as in an 



earlier paper [21] the tadpole character of the operator has not been recognized and 
its contribution to the parameter e' was strongly overestimated. The correct treatment of 



this problem has appeared in the most recent paper [fjj] which is in agreement with our 
analysis. 

Thus the problem is reduced to the estimation of the matrix elements of the operators 

(8) 

Q\ . At present there is no regular method to calculate that kind of object within 
QCD excepting the numerical simulations on the lattice. To estimate at least the scale 
of the non-leading l/m c corrections we will work with the simplified model. As a first 
approximation we take the operators that survive after the factorization procedure. Then 
one selects the operators containing scalar quark currents which can be written as 

{sLGf, u a^ u q R )(q R d L ) and {s L q R )(q R G ^v^di). 

This step seems to be justified because in the case of dimension six operators the 
similar "penguin-like" structures are strongly enhanced and dominate the others (see 
Sect. 3). The last simplification consists in the substitution 

-> m\qq (2.8) 

where mo determines the scale of non-locality of the quark condensate and is defined by 
the equation {qgsG^a^q) = ml(qq), ml(l GeV) = 0.8 ± 0.2 GeV 2 p3| . 



This substitution is valid in the chiral limit for the operator qg s G^ u a^ y q with dimension 
five in mass units. We suppose that it is justified also in our case at least for estimates up 

9 



to the order of magnitude. Actually all above assumptions as the factorization procedure 
in the case of dimension six operators become exact within the many color limit of QCD, 
N c — ► oo, to the leading order in iV c . 

Having adopted the assumptions described above the only operator which has the non- 

(8) 

zero matrix element for the considered process is the operator Q\ and the corresponding 
quantity reads 

(HQiV) = ^(HQ 6 |K>. (2-9) 

We should note that the coefficients C\ are finite to the leading order of the a s and 
independent of renormalization scheme. We can therefore use the leading order values of 
the mesonic matrix elements that is consistent up to the used level of accuracy. In the 
contrast, the next-to-leading a s corrections to Wilson coefficients depends on the renor- 
malization scheme and to make the physical amplitudes scheme independent, matching 
between Wilson coefficients and mesonic matrix elements in the same renormalization 
scheme has to be made [||. 

Thus, taking into account the first order l/m c corrections is reduced to the effective 
shift of the coefficients in front of the penguin operator Q$ in the effective Hamiltonian 

Using the numerical values z e = —0.015, ye = —0.102 at the point Aqcd = 0.3 GeV, 
jj, — 1 GeV , m t = 130 GeV |J one finds numerically the relative corrections to the 
Wilson coefficients in the form 

^-►^(l-O.l), y 6 ->y 6 (l + 0.01). 

The main correction appears in the real part of the Wilson coefficient of the penguin 
operator Q e . Parametrically, the contribution of the dimension eight operators can be 
as large as one half im^jvn^ ~ 0.5) of the one of dimension six operators, and not too 
small. In fact, we have found that when one estimates the mesonic matrix element 
of the local operator with dimension eight within the simplest factorization framework 
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the non- leading l/m c contributions to the kaon decay amplitudes are about 10% of the 
leading ones. However a large violation of the factorization for the matrix elements of the 
dimension eight operators does not seem to be impossible and the real value of the non- 
leading l/m c correction can be estimated only when somewhat self-consistent method to 
calculate these matrix elements within QCD will be available. 

3 The QCD penguin operator contribution to the 
K —> tttt decay amplitude. 



In this section we demonstrate a functioning of the method developed earlier in ref. [[L2 
for calculating a contribution of the penguin operator Qq to the K —>■ mi amplitude. This 
operator draws the special interest because its mesonic matrix element prevails over the 
current-current and other QCD penguin operators and gives a dominant contribution to 
the parameter e' (we do not consider the electroweak penguins). Moreover, the estimates 
of matrix elements of the current- current operators are much less controversial and are 
quite reliably given by the simple factorization procedure. We leave aside technical details 
of the calculations that can be found in ref. jl2| and focus on the main features of the 
approach. 

The matrix element of the K — > tttt decay amplitude due to the operator Q$ only 
reads 

(n + n-\Q 6 (fi)\K°) = (tt°tt°\QM\K ) = (tttt\Q^)\K°) . (3.1) 

For comparison with other approaches we use the following parametrization of this matrix 
element 

-(tttt\Q & (^)\K q ) = -B p f K m 2 K , f K = 1.23/ w , / w = 132 MeV (3.2) 

i 



where the Jk^k factor fixes a natural mass scale while B p is a dimensionless parameter 
to be computed. 

We start with constructing the chiral effective theory describing weak and strong 
interactions of pseudoscalar mesons at low energy H ^6f . The operator Q 6 is a composite 
quark operator which can be define accurately in the asymptotic freedom regime within 
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perturbation theory using the proper prescription of renormalization. Its matrix element 
can be hardly found in the practically interesting case of low energies where confinement 
takes place. However one can try to define the operator Q e at low energies within chiral 
perturbation theory, i.e. to find an effective realization of the operator Q e in terms of 
mesonic variables. The operator Q 6 belongs to the 8^ x 1 R irreducible representation of 
chiral SU(3) L x SU(3) R group and has the following SU(3) V flavor quantum numbers: 
S = 1, / = 1/2, ^3 = —1/2. To the lowest order of a chiral perturbation theory expansion 
there are two Lorentz invariant bosonic operators with relevant properties (d fM U ji d^U)23 
and (x^U + £/^) 23 where x stands for the pseudoscalar meson mass matrix and U = 
e -iV2<f>/U j g j-jjg un it ai y matrix describing the octet of pseudoscalar mesons. Hence the 
operator Q e can be represented as 

Q 6 = -ftWirfPV)* + 9'(x j U + U^ X )2z] (3.3) 

where g and g' are the dimensionless parameters which are not fixed by symmetry re- 
quirements. After such an ansatz the mesonic matrix elements of the operator Q e become 
exactly computable within the effective theory, in other words, we have found a suitable 
kinematical framework for our dynamical problem which consists in determination of g 
and g'. For our aim we need the matrix elements 

(nn\Q 6 \K°)=Mgf n m 2 K , (3.4a) 

(n + ( Pl )\Q 6 \K + (p 2 )) = f 2 w (4g(p lP2 )+4g'm 2 K ), (3.46) 

(0\Q 6 \K°) = -Aig'flm 2 K . (3.4c) 

We put the matrix element in eq. (3.4a) on the meson mass shell because this corresponds 
to the physical amplitude of K — > im transition. We also take into account that the 
term proportional to x is equal to a full derivative due to equations of motion and does 
not contribute to the amplitude at zero momentum transfer. The matrix elements in 
eqs. (3.4b) and (3.4c) are some auxiliary amplitudes considered only for determination of 
g and g' and can be computed at arbitrary point since they are momentum independent 
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to the leading order of the chiral expansion. So, we do not include the contribution of 
the so-called tadpole term (x U + U*x)23 to the physical amplitude (3.4a). As is well 
known the appearance of such a term is a consequence of working with a wrong vacuum 
solution. This term merely renormalizes the strong effective Lagrangian (beyond the mass 
shell) and can be absorbed into the meson mass matrix by a suitable SU(3)l x SU(3)r 
rotation. Thus, the tadpoles do not generate any observable effects for the effective chiral 
Lagrangian of order p 2 to the first order in Gf H 27 1. 



The problem is reduced to the determination of the single parameter g. It can be 
done by studying an appropriate three-point Green's function (GF) via the sum rules 
technique. We choose the GF in the form 



G„(p,q) =i 2 (U\Tjl{x)Q & {U)j\y)\uy^™dxdy 



= i 2 J (0\Tjl(x)Q 6 (y)f(0)\0)e^--^ + ^dxdy (3.5) 

where pi = p + q/2, P2 = p — q/2 and we take = d^^u, j 5 = uj^s as interpolating 
operators for pion and kaon fields. Because of dispersion relations, GF (3.5) looks like 
(j^ + (p2))(7r + (p2)\Qe\K + (p 1 ))(K+(p 1 )\f) R, R R 

2/2 2 \ ' 2 ' 2 2 ' ' ' ' \ / 

p 2 2 {pi-m 2 K ) p 2 , pi-m 2 K 

where the ellipsis stands for states without any kaon or pion poles, the quark currents 
projections on the mesonic states are 

(Ob> + (p)> = iUp» (K + ( P mO) = - i ^- (3.7) 

m s 

and we work with the massless pion (u-, <i-quarks). The representation (3.6) needs a 
comment. The matter is that to extract information about the matrix element in the first 
addendum in eq. (3.6) one has to distinguish the resonance contribution from miscella- 
neous ones. It is possible because, contrary to all other states, the resonance leads to the 
double pole in the dispersion relation. In representation (3.6) it can be achieved within 
the kinematics where (pq) = m 2 K /2. Setting (pq) = m 2 K /2 and multiplying eq. (3.6) by 
the expression (p 2 + q 2 /4 — m 2 K /2) one obtains 

(jf l \n + (p 2 ))(n^p 2 )\Q 6 \K+(p 1 ))(K+(p 1 )\f) 
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+ R L + R R + . . . (3i 



where the resonance contribution is explicitly distinguished as a pole term. Eq. (3.8) can 
now be treated by the standard sum rules technique. 

Straightforward calculations give the following asymptotic expansion for the GF 

G^(p,q) = PzpG(p, q) + . . . 

and 

Gfe *) = - (f) <« + ^ -> (f) <«+ 

+0(p~ 6 ) + (bilocal part) (3.9) 

where (uu) = (dd) = (ipip). We work in the leading approximation of the q 2 power 
expansion and simply put q 2 = 0. At the same time we have to keep the first order 
terms in the scalar product (pq) = m 2 K /2 and the strange quark mass m s and we also 
have to distinguish the strange quark vacuum condensate from the non-strange one 7 = 
(ss) / (ipip) -1/0 since this difference represents a parameter of SU(3)y symmetry 
breaking. As we see all terms of the zeroth order in m s cancel themselves in eq. (3.9) in 
agreement with the chiral structure of the operator Q§. The bilocal part of OPE does 
not contribute to the parameter g to the 0(p~ 6 ) order due to a specific non-symmetrical 
choice of the the interpolating currents in eq. (3.5) 

To determine the parameter g we apply the finite energy sum rules technique |28| to 
the function G(p 2 )(p 2 — m 2 K /2) which is defined by the equation 

G(p 2 ) = G(p,q) | 2 (3.10) 

q 2 =0, [pq)=-^- 

The result for the matrix element (tt + \Qq\K + ) reads 

(w mso)lK } = 2^^r s ° i 1 + i< J (3 - n) 

where Sq is a duality interval. Combining eq. (3.11) with the representation (3.4b) and 
using the PC AC relation —2m s {i(jifj) = f\vn? K we find 

5 + 2 9 = -T~2~r~n 1 + T~T ■ ( 3 - 12 ) 
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It can be easy shown |12| that the first addendum in brackets on the rhs of eq. (3.12) 
represents exactly the contribution of the constant g to the entire sum. Thus we have 

3 fx s 

9 = -^rr (313) 

The question now is which value for the duality interval Sq has to be used. Throughout 
the calculation we keep only the first order quantities with respect to SU(3)v symmetry 
breaking parameter - the strange quark mass; in eq. (3.11) the chiral suppression is 
present as a factor m s . Hence, for consistency, we should use for the duality interval Sq 
its value in the chiral limit, i.e. the value of the pion duality interval Sq = 0.8 GeV 2 
j29fl. There exists a prejudice that the duality interval in the pseudoscalar channel can 



be abnormally large due to a possible contribution of so-called "direct" instantons [30 



but the real estimate of that contribution is practically absent while the use of the above 
value reproduces the pion decay constant f n with reasonable accuracy. Actual value of 
the So parameter can be found only after adding the non-perturbative corrections due to 
higher dimension operators which are supposed to be small in our case. Finally for the 
parameter B p we obtain 

B P {sl) = ^r 2 % = 7- (3.14) 

J TV 

Eq. (3.14) is written for the normalization point /i 2 = Sq. For an arbitrary normalization 
point we have 

BM)=B P W(J^) (3.15) 

where 76 = — 14 is the anomalous dimension of the operator Q$. 

The strong dependence of the rhs of eq. (3.14) on the parameter sq is considerably 
smoothed by the renormalization group factor ( "[^°j ) 7//9 - Indeed, the result for the param- 
eter B p {l GeV) changes only from 8.19 to 9.38, i.e. less than 15%, when one substitutes 
the kaon duality interval Sq = 1.2 GeV 2 |2!| into eqs. (3.14, 3.15) (the QCD scale is 
chosen to be Aq C d = 0.3 GeV). 

Let us now estimate the uncertainty of our result. On the physical side the errors 
originate from the higher order terms in chiral expansion which are omitted in eq. (3.3). 
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Their relative weight can be represented by the ratio m 2 K /A 2 where A x is a chiral scale 



parameter. Both on theoretical and empirical grounds one expects [[31 



A 2 x = Sn 2 f 2 ~ 1 GeV 2 . (3.16) 

Thus, a choice of the representation for a quark-gluon operator (for example, Q e ) to the 
leading order in the chiral effective theory might introduce a 25% error in the estimate of 
the real physical matrix element. 

On the theoretical side of the sum rules the errors come from several sources: 

1. Corrections due to operators with higher dimensionality which start with the term 
(ipg 2 G 2 ip) (pq) I p 6 and, probably, some 0(p~ 6 ) bilocal contribution. They are under control 
and do not lead to any sizable variation of our results though their actual magnitude is 
questionable because it requires making some estimates of the vacuum expectation values 
of higher dimension and bilocal operator; 

2. Perturbation theory corrections to the coefficient functions of the leading operators. 
The correct inclusion of these correction requires a complete generalization of the effective 
Hamiltonian up to the a 2 order. Therefore we consequently work only with the first order 
in a s terms in the effective Hamiltonian and omit radiative correction to the non-leading 
operator Q§ since its coefficient function is already of the a s order. 

Thus, we have calculated the matrix element of the penguin operator via the QCD 
sum rules with the result 

-(7T7r|Q 6 (l GeV)\K°) = -(0.34 ±0.09) GeV 3 (3.17) 
i 

or, in terms of the parameter B p 

B p (l GeV) = 8.2 ± 2.1 (3.18a) 

where the error bars estimate contributions of higher orders of the chiral expansion. The 
analysis carried out above shows that the uncertainties coming from other sources are rel- 
atively small. This result can be important for the precise calculation of the e' parameter. 
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The value (3.18a) is somewhat larger (with a factor 1.3 when m s (l GeV) = 175 MeV) 
than the leading order result obtained within the 1/N C expansion framework |J 

However a sizable enhancement (factor 1.5 — 2.0) of the leading order result (eq. (3.18b)) 
due to the next-to- leading corrections in the 1/N C expansion has been found ||14j| . This 



observation is in good agreement with our calculations. Our result gets also into the 



interval given by lattice models [£L3 



B p {lGeV) = U±3 (3.18c) 



B p (l GeV) = 12 ( 17 **i e lX - ( 3 - 1M ) 



but it is somewhat smaller than optimistic estimate of ref. [J/J where the soft pion technique 
was used 

m s (l GeV) ) 

It has been also shown that the Wilson coefficient of the operator Qq is essentially 
increased by the next-to-leading a s terms @. For example, at the point Aqcd = 0.3 GeV, 
\i = 0.8 GeV, the numerical value of the leading order coefficient Zq'°'^ = —0.028 while 
inclusion of the next-to-leading a s corrections leads to the result z^" 1 '^ = —0.098. However 
the consistent generalization of the effective Hamiltonian beyond the leading logarithmic 
approximation requires to derive the corresponding mesonic matrix elements up to the 
same order in a s (see Sect. 2). Thus this enhancement in general can be canceled by 
radiative corrections to the matrix element of the operator Qq. If it does not occur 
the penguin operator becomes quite important in the analysis of AI = 1/2 rule since it 
provides approximately 20% of the physical value of the decay amplitude with the isotopic 
spin transfer AI = 1/2. 

4 The scalar gluonium contribution to the K —> tttt 
decay amplitude. 

As it has been already pointed out the possible source of the enhancement of the AI =1/2 
amplitude is the unfactorizable contributions to the matrix elements of the local four- 
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quark operators produced by the low-energy strong interaction of the light quarks. In 
the chiral Lagrangians approach these non-factorizable contributions reveal themselves in 
two different ways: first, they appear as corrections to the couplings characterizing the 
" factorizable" weak chiral Lagrangian in 0(p 2 ) and higher orders, second, some new non- 
factorizable terms emerge. The latter is the case for K — > tttt decay mode with gluons 
forming an intermediate state. 

The unfactorizable corrections of the first and of the second type have already been 
discussed in the framework of 1/N C expansion [|| ||. In ref. ||] the non-perturbative 
corrections to the weak chiral Lagrangian of the leading 0(p 2 ) order caused by the presence 
of the non-perturbative gluon condensate was derived while in ref. || the 0{j> A ) order 
contribution of the unfactorizable chiral loops was computed and the sizable violation of 
the factorization has been found. However in both pointed approaches the above mode 
was omitted. Meanwhile it is quite reasonable to suppose that the gluons being considered 
as exchange particles can be important in explanation of the AI =1/2 rule because they 
carry a zero isospin and contribute only to the AI = 1/2 amplitude. Moreover that kind 
of contribution being non-leading in chiral expansion nevertheless can be sizable due to 
the strong effects of the non-perturbative gluon vacuum. 

Before treating the long-distance effects of the meson-gluon transitions it is very useful 
to consider the similar phenomenon arising already in perturbative QCD as a leading cor- 
rection in the inverse mass of the charmed quark given by eq. (2.3). If we restrict the anal- 
og 1 ) r*7\ 

ysis to the scalar colorless gluon configuration G a ^ v G a ^ v only the operators Q 3 , m s Qi , 
and m s Q^ give a contribution and eq. (2.3) gets the form 

HP = ^V ud V: s (l - r) (-^™ S ^7^G;j + h.c. (4.1) 

To derive the effective chiral Lagrangian involving Goldstone degrees of freedom only 
which then can be used for the calculation of the decay amplitude one has to replace the 
QCD operator 

m s s R d L -^Gl u Gl v (4.2) 
in eq. (4.1) by its mesonic realization. Using its chiral transformation properties one can 
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write down the representation 

m s s R d L ^G^G% = Aft(U\) 2 , + Bf^xh3tr fl (d^d»U) + 

+ (other 0{p A ) terms) + 0(p 6 ) (4.3) 

where A and B are the dimensionless parameters. The second term on the rhs of eq. (4.3) is 
separated from the other 0(p 4 ) structures for the reason will be clarified below. The 0(p 2 ) 
term in eq. (4.3) is exactly a tadpole term which should be omitted. So the contribution to 
the physical amplitude is determined by the 0(p 4 ) part of the chiral representation (4.3). 
The most transparent way to obtain this part is to consider the quark-gluon operator 
(4.2) as a product of the (pseudo)scalar quark current and scalar colorless gluon operator. 
Then one can replace the quark current by its mesonic realization according to the PC AC 
hypothesis and current algebra [26[] 

m 9 s R d L ^-S.{U^ X )23. (4.4) 

o 

On the other hand there is a low-energy theorem based on the fundamental properties of 
the energy-momentum tensor which gives the chiral representation of the gluon operator 



32, 33 



^G%G«„ = -jjltr fl {d^d»U) + 0{p% (4.5) 



Eqs. (4.4, 4.5) allow us to determine the parameter B 

which now modulates the unique term of 0(p 4 ) order in eq. (4.3). This approximation 
corresponds to the simplest physical picture where the kaon is annihilated by the pseu- 
doscalar quark current while the pion pair is born by the gluon operator. We should not 
that the above separation of the operator (4.2) into the quark and gluon parts leads to 
some uncertainties in the final result since the effects of the interaction between them 
are lost. However the corresponding corrections seem to be suppressed at least at the 
perturbative level. Eqs. (4.1, 4.3-4.6) result in the effective chiral Lagrangian of the form 

L c G h = -^V ud V: s (l - r) [-^^ X ) 2 ,tr fl (d^d^U)] + h.c. (4.7) 
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Since the pion pair is born by the gluon operator this Lagrangian describes the investigated 
decay channel with gluons forming an intermediate state. 

Now the corresponding K — > rnr decay amplitude becomes explicitly calculable. We 
will use the standard parametrization of the amplitude A with the isospin transfer AI = 
1/2 

ReA = —j= sin 6 C cos 9 c g 1/2 f K m 2 K (4.8) 

where 9 C stands for Cabibbo angle, gy 2 is a dimensionless parameter. The new contribu- 
tion reads 

I Tfl 2 

A 9i/2 = ^r-f ~ ID" 3 . (4.9) 
30p 3 m l c 



While the experiment gives |L9| 



^ = 3.9 (4.10) 

and the most recent theoretical estimation is gy 2 ~ 2.4 0]. As we see the local (pertur- 
bative) part of the new decay mode is negligible according to the general estimate of the 
scale of the leading order charmed quark mass corrections done in Sect. 2. 

However the local effective Hamiltonian (4.1) does not exhaust the whole physics of the 
meson-gluon transitions. It cannot account for the long-distance contribution connected 
with the propagation of the soft w-quark round the loop of the annihilation diagram. 
Because of the lightness of the w-quark this contribution can not be represented as a 
local vertex. It ultimately depends on the infrared properties of QCD and requires non- 
perturbative approach. We will follow the general line of the approach applied in Sect. 3 
to the calculation of the penguin operator matrix element. 

In so doing we start with a tree level Hamiltonian which after decoupling of the c-quark 
has the form 

HXs=i = %V ud V: s Q 2 + h.c. (4.11) 

The quantity of interest is an effective theory realization in terms of mesonic variables of 
the part of the operator Q 2 = Q 2 which is responsible for the kaon transfer into gluons. 
Invoking the results of our previous consideration we can write down this part in the form 

0,2 = g G fl{tf X )2ztr fl {d^d»U) (4.12) 
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where g is a dimensionless parameter to be computed. We should note that the chiral 
representation of the whole operator Q2 contains a large number of structures but we 
are interested only in the part corresponding to the transition with the gluons forming 
an intermediate state which has the unique representation (4.12). Indeed, there is the 
single SUy(3) octet terms (4.12) in the chiral weak Lagrangian describing the K — > tttt 
decays which is proportional to the trf^d^U^d^U) as it required by eq. (4.5) ||. Thus 
the problem is reduced to the computation of the chiral coupling constant g G . As for the 
parameters g and g' of the chiral representation of the penguin operator it can be done 
by studying the appropriate GF via QCD sum rules technique. For the technical reason 
working with a two point GF is preferable. In the given decay channel the pions are born 
by a gluon cloud therefore the gluon operator G a G a can play the role of an interpolating 
operator of the pion pair. Thus, it is natural to choose GF in the form 

G(p) = J (0\T^G%G%(x)Q 2 (0)\K°(q))e^dx\ q=0 . (4.13) 

A remark about the chiral limit for the kaon in eq. (4.13) is necessary. The repre- 
sentation (4.12) fixes the correct 0(p 4 ) chiral behavior of the considered decay amplitude 
and does not depend explicitly on the kaon momentum. Keeping a non-vanishing kaon 
momentum leads to a shift of the decay amplitude that lies beyond the accuracy of the 
present approach. Thus we can put q = in eq. (4.13) and work with GF depending on 
one argument only. 

Saturating GF (4.13) by the ir + n~ and states (the lowest states with the proper 
quantum numbers), substituting the operator Q2 for its mesonic realization (4.12) and us- 
ing the low-energy theorem (4.5) one obtains at small momentum p the following physical 
representation 

Gip) = 9 G ^r^p 4 \n (^) + 0(p% (4.14) 
The theoretical side reads after making use of OPE 

G(p) = i-L l n ^(0\m s s R g s GlX^d L \K°(q))\ q=0 + 0(a 2 p 2 ) + 0(p- 2 ). (4.15) 
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Factor m s in eq. (4.15) provides the correct chiral property of GF and justified the repre- 
sentation (4.12) for the operator Q 2 - By contraction of the kaon state one can transform 
eq. (4.15) into the expression 

G(p) = ^[=f)^f K rnWo- (4-16) 

For extracting information about the chiral coupling constant g G we use finite energy 
sum rules with the result 

G _ Zfofic flml a s (s ) 

128 J n S IT 

To take into account the strong interaction at short distances the operator in the 
effective Hamiltonian has to be multiplied by the corresponding Wilson coefficient ^(so)- 
Finally the new contribution to the theoretical estimate of the A J = 1/2 amplitude in 
terms of parameter gi/ 2 takes the form 

i\gi/ 2 - z 2 (s ) — -J — — • (4.18; 

O 6q 7T 

This result needs some comments. 

1. This next-to-leading in 1/N C expansion contribution is missed within the factorization 
framework and also within any approach where quark currents in four-quark operators 
are replaced by their mesonic counterparts separately. 

2. The gluon cloud in the intermediate state does not form a resonance state and, there- 
fore, the contribution (4.18) is not suppressed by a large scalar meson mass. 

3. In general the more complicated scalar colorless gluon configuration, for example 
f abc G a ^ v G b vX G c XfjL) could be an intermediate state in this channel as well. However the 
theorem (4.5) shows that the two pion form factor of these configurations can be of 0(p 4 ) 
or higher order in chiral expansion and that lead only to the negligible 0(p 6 ) shift of the 
decay amplitude. Indeed, the theorem (4.5) sets the equivalence between the trace of the 
energy momentum tensor expressed in terms of QCD and mesonic degrees of freedom. 
In chiral limit in the leading order in momentum expansion there is a single Lorentz 
invariant flavor singlet mesonic configuration which is proportional to the trace of the 



22 



energy momentum tensor. On the other hand there are no extra gluon contributions to 
the conform anomaly. 

The question now is what numerical value for the duality interval so has to be used. 
Actually, the allowed value of the duality threshold is quite restricted by the form of the 
physical spectrum and by the requirement of absence of uncontrollable a s corrections. To 
suppress contributions of higher mass states, for example, a scalar meson cr(0.9 GeV), to 
the considered channel one has to take sq < (0.9 GeV) 2 . At the same time the physical 
representation (4.14) is obtained in the leading order in chiral perturbation theory and the 
whole procedure is justified until the ratio So/A^. remains small. On the other hand at the 
scale \i less than 0.8 GeV the perturbative a s corrections to Wilson coefficients become 
uncontrollable [0 and for consistency of the approach one has to set the low limit of the 
duality interval to be so > (0.8 GeV) 2 . The reasonable choice for the duality interval now 
reads s = (0.8 GeV) 2 . 

Following the approach of the Sect. 3 let us estimate the uncertainty of our result. On 
the physical side of sum rules the errors related to higher order terms in chiral expansion, 
which have been omitted in eq. (4.12), are, in general, unknown. But one can hope that 
in the spirit of chiral perturbation theory they are about 25% J2(| |^T|. On the theoretical 
side of the sum rules the errors come from two sources. The first one is the perturbative 
part of OPE (the unit operator) that is suppressed by a loop factor a s j^K ~ 10~ 3 and 
cannot lead to a sizable change of our result. Next non-perturbative corrections due to 
operators with higher dimensionality seem to be more important. They start with the 
dimension eight operators which have already been discussed. Numerical estimates would 
require knowing the matrix elements of those operators between the kaon and the vacuum 
state which are not available now. But as a first approximation the relative weight of these 
corrections can be represented by the ratio 



T(5) ' 4vr 2 r(3) 



0.1 (4.19) 



where {^G a G a ) ~ (330 MeV) [M i the Gamma function factor 1/T(n) comes from 



the quark loop with the n — 1 gluon field or mass insertion. Here we find the same 
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situation as in the analysis of l/m c corrections where the contribution of the dimension 
eight operators is suppressed rather numerically than parametrically. 

Taking into account the uncertainty of determination of the parameter m§ we estimate 
the error bound to be about 40%. Numerically one obtains 



at the point s = (i 2 = (0.8 GeV) 2 , A QCD = 300 MeV, z 2 (s ) = 1.49 [|. 

Thus, the new contribution provides about 15% of the experimentally observable am- 
plitude (4.10). At the same time it is comparable with the leading order result for the 
decay amplitude obtained by the naive factorization of the four-quark operator Q 2 when 
all strong interaction corrections are neglected 



This implies the strong violation of factorization in the 0(p A ) order in chiral expansion 
that leads to additional enhancement of the theoretical estimate of the K — > tttt decay 
amplitude with the isospin transfer AI = 1/2. 

5 Conclusion. 

In the present paper we make an attempt to clarify the role of annihilation or "penguin" 
mode in the descriptions of the K — > tttt decay in the Standard Model. We concentrate 
our attention on some new aspects which have not been considered yet. 

The complete form of the effective low-energy AS = 1 Hamiltonian up to the first 
order in m s , l/m c and a s with m t ~ M\y is obtained. The violation of factorization for 
the mesonic matrix element of the operator Q 2 caused by K — > tttt transition with the 
gluons in the intermediate state is also analyzed. 

We propose the regular method to evaluate the mesonic matrix elements of the effective 
Hamiltonian for K — > tttt transitions. The weak Hamiltonian is represented through 
the variables of chiral effective theory and the coefficients of proportionality between 
quark-gluon operators and mesonic operators are determined via QCD sum rules for an 



Ag 1/2 = 0.56 ±0.22 



(4.20) 



9[J 2 = 5/9. 



(4.21) 
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appropriate three-point GF. We apply this technique for computing the QCD penguin 
operator matrix element. 
The main results are 

1. The contribution of the dimension eight operators (non-leading l/m c corrections) 
proves to be small numerically (about 10 % of the leading term) when one estimates 
the corresponding mesonic matrix element within the simplest factorization framework. 
However its real magnitude can be estimated only when more detailed information on 
the matrix elements of the local operators with dimension eight becomes available, for 
example, from lattice calculations. 

2. The quark-gluon operator m s SRg s G^ u a^ u d L contributes to the imaginary part of the 
effective Hamiltonian in the case of a heavy top quark already in the first order in a s due 
to incomplete GIM cancelation. However its contribution to the parameter e' is suppressed 
because of the smallness of the corresponding Wilson coefficient and the specific (tadpole) 
chiral structure of this operator. 

3. Sum rules predict for the matrix element of the QCD penguin operator 

-<7T7r|Q 6 (l GeV)\K°) = -(0.34 ±0.09) GeV 3 
i 

where the error bars estimate contributions of higher orders of the chiral expansion. The 
analysis carried out above shows that the uncertainties coming from other sources are 
relatively small. This value is somewhat larger than the results obtained within the 1/N C 
approximation framework j|, || but it is somewhat smaller than optimistic estimate of 



ref. 0. Our result gets also into the interval given by lattice models [fT3~ |. The relatively 
large value of the matrix element along with the strong enhancement of the corresponding 
Wilson coefficient [[| make the penguin operator quite important in the analysis of the 
AI = 1/2 rule since it provides about 20% of the decay amplitude with AI = 1/2. 
4. The non-perturbative contribution to the mesonic matrix element of the leading op- 
erator Q2 produced by the transitions with gluons in the intermediate state gives an 
additional enhancement of the theoretically predicted K — ► tttt decay amplitude with 
AI = 1/2 and provides about 15% of the experimentally observable amplitude value. New 
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contribution is of 0(p 4 ) and is lost within the factorization framework but numerically 
is comparable with the result of naive factorization for the AI = 1/2 decay amplitude. 
It allows us to conclude that there is a sizable violation of the factorization in the 0(p 4 ) 
order in chiral expansion. 

To conclude, the comparison between the "penguin-like" and "non-penguin" parts 
of the AI = 1/2 amplitude is quite instructive. According to the recent work ]2[] the 
current- current operators Q\ and Q2 (without Q 2 part) give about 45% of the experi- 
mentally observable amplitude while the total effect of the "penguin- like" contributions 
of the operators Q 6 and Q 2 turns out to be about 30 — 40%. So large "penguin-like" 
contribution improves considerably the theoretical description of the AI = 1/2 rule in 
non-leptonic kaon decays. Indeed, our results imply ReAff ~ 0.8ReAQ Xp while ref. 
gives ReAff ~ §.QReA £ Q P with penguin contribution being only about 15%. As we see, 
the discrepancy between theory and experiment still exists but it is getting smaller. This 
result can be presented in more impressive form if one considers the ratio ReA /ReA 2 
where A2 stands for the amplitude with the isospin transfer AI = 3/2 [^]. Then we obtain 
{ReA /ReA 2 ) th ~ 21. that is very close to the experimental data (ReA /ReA 2 ) exp = 22.3. 
We should emphasize that this result is obtained without the standard practice of using the 
extremely low normalization point or extremely light strange quark. Since the resources 
of the perturbation theory seem to be exhausted the further progress in the theoretical 
explanation of the AI =1/2 problem will probably be connected with studying the un- 
factorizable contributions to the mesonic matrix elements caused by strong interactions 
at low energy. 
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